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A model of superconductivity in layered high-temperature superconducting cuprates is proposed, based on 
the extended saddle point singularities in the electron spectrum, weak screening of the Coulomb interaction 
and phonon-mediated interaction between electrons plus a small short -range repulsion of Hund's, or spin- 
fluctuation, origin. This permits to explain the large values of Tc, features of the isotope effect on oxygen and 
copper, the existence of two types of the order parameter, the peak in the inelastic neutron scattering, the positive 
curvature of the upper critical field, as function of temperature etc. 
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PACS numbers:74.20.Fg,74.25.Kc,74.62.Dh,74.72.Bk 
I. PHONON-MEDIATED INTERACTION 

The origin and the symmetry of the order parameter in the 
HTSC is the primary topic in the physics of high-temperature 
layered cuprates. The crystalline symmetry, even if it is as- 
sumed tetragonal, and the order parameter to be real (this 
follows from time-reversal symmetry), requires only that at 
a 7r/2 rotation in the ab-plane it either does not change, or 
changes its sign. Symmetry by itself does not tell us how 
many nodes must exist in either case, or whether the for- 
mula cos kx — cos ky is correct. Therefore the only way 
to find the order parameter is to study the interaction be- 
tween electrons, which leads to pairing. The high values 
of Tc and absence of the isotope effect in optimally doped 
Y Ba2Cu^O'j-s (YBCO) led to the conclusion that this in- 
teraction cannot be mediated by phonons. However the ex- 
periments by the groups of J. Franck and D. Morris [jl|,^ 
showed that a partial substitution of Y by Pr, or of Ba 
by La leads to the increase of the isotope effect simultane- 
ously with the decrease of the critical temperature, so that at 
Tc « 40X0: w 0.4 (Tc cx M"") (Fig.l). An alternative 
appears: either these substituted compounds have a differ- 
ent mechanism of superconductivity compared to YBCO with 
Tc ~ 90K, or the mechanism is always phonons, and the 
absence of the isotope effect in YBCO is due to something 
else. I prefer the last option. In its favor speak also the tunnel- 
ing spectra at higher voltages in Nd2-xCexCuOi (NCCO) 
g (Fig.2) and in Bi2Sr2CaCu20s (BSCCO) [|] (Fig.3). 
The Eliashberg interaction constant a^Fiuj) extracted from 
the tunneling conductance has a shape corresponding to the 
phonon density of states F{uj) obtained from inelastic neutron 
scattering. The deep minima in a^F(ix>) compared to F{uj) 
can be attributed to the fact that electrons interact differently 
with different phonons. Lately the strong electron-phonon in- 
teraction in BSCCO was confirmed by appearance of branch- 
ing on the current-voltage characteristic of a break junction at 
eV < 2A [^], representing phonon-assisted Josephson effect 
(similar to Shapiro steps) for low-frequency optical phonons. 
The frequencies fitted the Raman data (see sec. 6). 



II. EXTENDED SADDLE POINT SINGULARITIES 

A few years ago two experimental groups: J.-C. Cam- 
puzano's from Argonne ||^] and Z.-X.Shen's from Stan- 
ford [^, investigating angle resolved photoemission spec- 
tra (ARPES), discovered that the dependence of the energy 
of quasiparticles on their momentum e{kx,ky) (the depen- 
dence on kz is weak due to quasi-two-dimensionality) has 
flat regions, where it does not depend on one of the compo- 
nents, so that the spectrum becomes quasi-one-dimensional 
in these regions (Fig. 4). This was first discovered in YBCO, 
YBa2Cu40s, BSCCO, and appeared also in band structure 
calculations by several groups, particularly, A. Freeman's ||^], 
for mercury compounds. Such a feature in the spectrum was 
called "extended saddle point singularity". The density of 
states in the vicinity of such a singularity is proportional to 
{e-€n)-^/^. If the Fermi energy is close to eg, this enhance- 
ment is more than sufficient to explain the high values of Tc 
assuming the usual strength of the electron-phonon interac- 
tion. 

The second feature is that the integral entering the BCS 
equation 



A(p)= J K{p,p')f[A{p') 



(27r)3 



(1) 



which is logarithmically divergent, if the density of states is 
constant, becomes convergent. In the case, if ep — eo <C 
ujQ, ojo being some characteristic phonon frequency, the criti- 
cal temperature does not depend on ujo, and hence there is no 
isotope effect. The BCS-type formula for Tc for this case has 
the form: 



87(eF - Co) 
Tc = exp 



1/2 2 

gnix K, 



(2) 



The most important here is the appearance of e^? — eo in the 
exponent and in the prefactor (I will not explain the meaning 
of other notations in this expression). If the distance ep ~ 
eq increases, Tc decreases, and, at the same time, cf — eo 
in the prefactor is gradually replaced by wq , enhancing the 
isotope effect [P] . This fits the experimental data. I would like 
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to stress that the usual saddle point singularity, extensively 
studied by several groups, e.g. J. Friedel, A. Labbe, J. Bok, 
D. Newns, R. Markievicz (see | ]lO| ] for references) does not 
lead to such consequences. The enhancement of Tc is weak, 
and no new energy scale is introduced in the BCS integral, 
so that the absence of the isotope effect cannot be explained. 
Therefore the statement that the flat regions are equivalent to 
2D van Hove singularities is wrong. 



III. BAD COULOMB SCREENING 

Under proper assumptions the flat region can dominate in 
the total density of states; then, however, we encounter an- 
other puzzle. The ARPES experiments of the same groups es- 
tablished a marked anisotropy of the energy gap (in the BCS 
theory the gap corresponds to the absolute value of the or- 
der parameter). The usual belief is that with phonon-mediated 
interaction no anisotropy of the order parameter can be ob- 
tained, even in case if the energy spectrum of the metal is 
anisotropic. This is due to exchange of phonons with high 
momenta which effectively connect all points of the Fermi 
surface in the BCS equation (1) (see Fig. 5). However, this 
statement, being true for good metals, is not so evident for 
the substances under consideration. According to numerous 
expe-rimental data, they are intermediate between metals and 
ionic crystals. The cross-over between these two limiting 
cases can be imagined as the change of the Debye screening 
radius from interatomic distances to infinity, and the HTSC 
can be viewed as substances where all Coulomb-based inter- 
actions are screened at large distances. This refers also to the 
electron attraction mediated by phonons, since it is due to the 
electron-ion interaction. Therefore it is reasonable to assume 
this interaction in the form 



V{k) 



(e-e')'-C.2(fc) 



(3) 



where k K is the reciprocal Debye radius, ^ is the elec- 
tronic energy, and K - the reciprocal lattice period. Since 
different models lead to different powers n, we will simply 
assume n > 1. Under such circumstances the exchange of 
phonons with small momenta becomes dominant, and we can- 
not reach the singular region from some distant place (Fig. 5) 
without decreasing drastically the value of the order parame- 
ter. Therefore it becomes very anisotropic, and its maxima are 
located in the regions with the maximal density of states JTl]]; 
this, in fact, is observed in experiment. Moreover the values, 
distant from the singular regions, will be defined by the max- 
imal value, the connection being established through the BCS 
equation. Therefore the critical temperature is defined only by 
these singular regions. 



IV. TWO TYPES OF THE ORDER PARAMETER 

In a purely phonon picture there can be no changes of sign 
of the order parameter, and this contradicts many reliable ex- 
perimental data. However, apart from phonons, there can exist 
interactions of some other origin. The isotope effect tells us 
that the phonons dominate. Therefore these other interactions 
can become important only, if their nature differs substantially 
from the phonon-mediated attraction. We will assume that this 
additional interaction is a short-range repulsion coming either 
from spin fluctuations, or being the renormalized repulsion at 
the Cu sites. To the first approximation different regions with 
the maximal density of states and the maximal absolute value 
of the order parameter are disconnected. Therefore, the order 
parameter can be either even, or odd with respect to a tt/2 
rotation in the ab-plane. The corresponding formulas are: 



A(^) 2(n-l)Pc 



Oy 




Ai 



A{ip) _ 2(n - l)Poy ( n 



U 
9 _ 

2n 



(4) 



Ai 



sm 



■(I) 



(5) 



where Ai is the order parameter in the singular region. 
Pay is the length of the singularity, d is the period in the c- 
direction, po is the radius of the Fermi surface, which we as- 
sume cylindrical, and ipo = tt/2 are the locations of singu- 
larities, and U is the amplitude of the short-range repulsion, 
which we suppose isotropic. These expressions are valid far 
from the singular regions. A simple interpolation permits to 
describe the whole Fermi surface. The simplest formulas for 
n = 1 are 



A, 



A 



Ai 

^odd 



A+1- cos(4(/?) ' 

ylcos(2iy9) 
' ^ + 1 - cos(4(/7) 



(6) 
(7) 



(A-adjustable constant), representing the even and odd cases. 
Fig. 6 presents the plot of eq. (^ for A = 1. The order pa- 
rameters (4) have the same symmetry, as the"s-wave" and "d- 
wave" although they differ from the usual expressions. As al- 
ready said, to the first approximation the even and odd states 
are degeneratre. This degeneracy is lifted, if we take into ac- 
count the energy of the small overlap of the order parameter 
originating from different singular regions. In case, if there 
were only phonon attraction, the even solution would have 
definitely the lowest energy. However, the presence of a small 
short-range repulsion with an isotropic Fourier amplitude in 
momentum space can drastically change the situation JT^]. 
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The interaction Hamikonian in momentum representation has 
the shape: 



PiP2k 

^(Pi - P2)a/3(-p2 + k)aQ(p2) 



U 



+ 2" ai(Pi)al3(-Pi +k)a/3(-P2 +k)aa(p2) (8) 

PiPzk 

where F(k) — [^^/(fc^ + k^)]". Since the main contribu- 
tion comes from vicinities of the singular regions, which are 
in the a- and b- directions, we separate parts, associated with 
the condensate by putting fc = and transforming 

(9) 

where lap is the unit antisymmetric matrix: / ~ iay. The 
can be considered as a creation operator of a Cooper pair, 
consisting of electrons with momenta pi and — pi. We leave 
quadratic and linear terms in F^. Then we get 



Hi, 



2P, 



yO 



{2T:yvid 



[Al($a + *h)+2A2$a$b] 



Ai 



*a/a/3(a/3(p)aa(-p) + al{-p)a}j{p)) + 



pGa 



<^>blaf}{apip)aj-p) + al{~p)al{p))] 



peb 



-[^ ^ala0iapip)aai-p) + flj, (-p)a|3 (p) ) 



peb 



J2 *b^"/3(a/3(p)aa(-p) + al{-p)al{p))] (10) 



where 



pGa 



Ai 



Ai = 




(11) 



(12) 



We used here that in the singular region the functions Fa,b de- 
pend only on the normal component of the momentum, i.e., on 
£, — vi {pn — Po), and performed integrations over other com- 
ponents. Since we assume that gn^ 3> UPyo/d and k <C Po, 
we have: Ai ^ A2, and hence the overlap terms in the Hamil- 
tonian are indeed small. We would like to stress here that the 
reduced interaction Hamiltonian (p^, where only the singu- 
lar momentum regions were left, is sufficient for calculation 
of the critical temperature and total energy but fails when the 



low-energy excitations play the major role, as in the specific 
heat and kinetics at low temperatures. Applying all the super- 
conducting machinery (see ref. [|l3|]) to the full Hamiltonian, 
including kinetic terms, we obtain the necessary quantities: 

T'c.odd = 8-^1 exp[-l/(Ai - A2)] 



7 



Aiexp{-A2/[Ai(Ai -A2)]} 



7 



T 

c.o 



dd 



Tr_. 



« ^Aiexp(-A2/A?) 

TT 

Aiexp(A2/A?), A2 < Ai 
« -2^Aisinh(A2/A?), 

TT 

Ai =8Aiiexp(-l/Ai) 



(13) 



f^odd(o) - n 

even (0) 



4P, 



yO 



(27r)2uid 



A2sinh(A2/A2) (14) 



From (H) and it can be seen that, if A2 > 0, the equilib- 
rium solution is even, and, if A2 < , it is odd. This occurs 
at all temperatures, where superconductivity exists. From the 
definitions ( [l3| ) it follows that the phonons dominate in the 
interaction, if 

gKy{n~l):^2UPyo/d (15) 
Under this condition A2 can have either sign. It is negative, if 



2MJi 



U > 5[«7(2Po)] 



(16) 



The main conclusion is the division of layered cuprates in 
two very distinct groups. Strong anisotropy, originating from 
extended saddle point singularities and a long screening ra- 
dius, makes it plausible that the order parameter becomes odd, 
whereas weak anisotropy leads to an even order parameter 
Since the quasi- ID energy spectrum leads to a singularity in 
the density of states only from one side, it is likely to be ab- 
sent in the n-type cuprates, as NCCO . This can be the origin 
of their lower Tc and of the even order parameter 



V. SCATTERING FROM IMPURITIES 

A strong proof of the odd order parameter is the action of 
impurities [p^. Regular substitutional impurities are most 
likely ionized, and being small-angle scatterers, they do not 
influence Tc substantially. Of course, there is also some large- 
angle scattering but it has a small probability. The usual esti- 
mates lead to the conclusion that the critical temperature is in- 
fluenced only if rep ^ 1, which could be considered in favor 
of s-pairing. The fault of such an estimate is that a much larger 
T should be taken, namely, that of large angle scattering, and 
it could fit the condition tA '-^ 1. In the experiment [p3|], de- 
fects were created in YBCO by electron bombardment. The 
electron energy was chosen so that the oxygen atoms were 
displaced but not knocked out; this produced dipolar defects 
causing large-angle scattering. They were very effective in 
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suppression of the critical temperature. At the same time the 
copper atoms were not involved, and this excluded magnetic 
defects. The observed dependence of Tc on the concentration 
of these defects fits the prediction for an odd order parameter: 

where is the scattering probability proportional to the 
residual resistance in the normal state, and 

- (^/47)ri^'-") (18) 

is the critical value, where Tco vanishes. In reality, however, 
this does not happen, because the metal can undergo a phase 
transition into the even phase. This was not observed in the ex- 
periment [|l5]| in the accessible temperature interval, however 
there exist other ways to suppress the d-wave order parame- 
ter, and at least in one of such experiments performed by the 
group of L. Greene a "subdominant" s-wave order parameter 
was indeed observed [Q. 

VI. NEGATIVE ISOTOPE EFFECT ON COPPER 

Rather close considerations can explain the negative iso- 
tope effect on copper, discovered by the Franck's group Jl7|]. 
It was found that in underdoped samples of YBCO the substi- 
tution ^^Cu — Cu causes a change of Tc corresponding to a 
negative isotope effect, namely, a larger copper isotope mass 
leads to an increase of the critical temperature. There exist 
some other examples of a negative isotope effect, e. g., re- 
sulting from a substitution of hydrogen by deuterium in solid 
solutions. The distinction of the present case is that the same 
substances have a positive isotope effect for the oxygen sub- 
stitution ^^O — ^^O. The absolute value of the negative copper 
isotope effect depends on the oxygen concentration. In opti- 
mally doped samples it is small, otherwise it depends non- 
monotonically on concentration. The maximal absolute value 
is around 0.4; we have in mind the power a in the dependence 
Tc cx M~", or for small variations of M, 

a = -A(hiTc)/A(hiAf). (19) 

In a simple model one cannot hope to reproduce all the de- 
tails, and our goal will be only the qualitative explanation of 
the most spectacular features. All phonons, whatever their na- 
ture is, contribute to the attractive interaction, and therefore, in 
order to explain the negative isotope effect, we have to think 
about some different action. In this connection we may re- 
member that in a d-wave superconductor ordinary potential 
nonmagnetic scatterers with large angle scattering suppress 
superconductivity. Let us suppose that there exist phonons 
which do not affect strongly the overall ionic charge. The in- 
teraction of electrons with these phonons will be short-ranged, 
the same as with neutral defects. If we suppose also that the 
characteristic frequency of such phonons is small compared 



to Tc, we obtain a situation similar to elastic impurity scatter- 
ing. In a d-wave superconductor this scattering will suppress 
superconductivity, and its effect will increase with increasing 
amplitude of oscillations, i. e. with decreasing mass of the 
ions involved. This consideration is indeed justified by exact 
calculations. The result is 




where 1^02 and are the frequency of Cu phonons with 
momentum po\/2 and their interaction constant with elec- 
trons, V is the density of states. The effect increases in magni- 
tude with decreasing Tc, i. e. with underdoping. For an s-type 
order parameter the isotope effect is always positive. The Ra- 
man scattering experiments |jT^], [|l9p show that indeed low 
frequency phonons do exist. The best candidate is the mode 
with cjQ w 150cm^^ which corresponds to vibrations of Cu 
atoms. For this mode the expansion parameter of the theory 
(wo/27rTc)2 - 0.24 (we took T^ = IQK). 

VII. NEUTRON PEAK 

Another interesting consequence of the proposed theory is 
the explanation of the origin and the prediction of the fine 
structure of the peak in inelastic neutron scattering in YBCO 
||2C|]. The clear narrow peak was observed in the super- 
conducting phase for momentum change close to (7r/2, 7r/2) 
and energy change close to 41 meV ( [pT|], and references 
thereon), which corresponds, approximately, to 2A,„a^, as 
obtained from tunneling measurements p^ ] . The momentum 
change corresponds to the transition from one extended saddle 
point singularity to another one. In these regions the spectrum 
is almost one dimensional. At Fig. 7 the energy in the normal 
state is presented as function of the corresponding momen- 
tum: for the a-vicinity, and py for the b-vicinity (dashed 
curve). In order to obtain from here the superconducting spec- 
trum one has to reflect the curve with respect to the chemical 
potential (straight horizontal line) and introduce gaps at the 
intersections. Then we get the two continuos curves at Fig. 7. 
The momentum and energy lost by the neutron are absorbed 
by the electron system and can be described, as transitions 
from the lower to the upper curve. Singular points in the scat- 
tering probability as function of energy can appear, if both, 
the initial and the final point correspond to the maxima of the 
density of states, i.e. to the extrema of the curves. The lowest 
energy corresponds to the transition between the maximum of 
the lower curve (1) to the minimum of the upper (1'), and at 
T = this is the threshold, located at 2D. Then comes the 
transition from the maximum of the lower curve (1) to the 
maximum of the upper (2'). Ths must be a maximum but it is 
not a threshold, since processes with smaller and larger ener- 
gies are possible. Finally comes the transition from the mini- 
mum of the lower curve (2) to the maximum of the upper (A'). 
This corresponds to the end of an absorption path, and hence 
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the scattering probability must have a drop at this point. Rig- 
orous calculations confirm this picture. The result is presented 
at Fig. 8. This has to be compared with the experimental curve 
at Fig. 9 [pO|]. It should be added that due to the so called 
"coherence factor" the maximum appears only, if the order 
parameter has different signs for the initial and final states (d- 
wave symmetry). This makes the inelastic neutron scattering, 
as well as the negative isotope effect described above rather 
unique bulk measurements confirming the d-wave order pa- 
rameter. The measurements associated with surfaces and in- 
terfaces can reveal properties different from the bulk, as, e.g., 
the suppression of the d-wave order parameter and appearance 
of a subdominant s-wave. Another important feature is that in 
a two-layered substance, which is YBCO, the maximum ap- 
pears only in the channel which corresponds to the interaction 
of neutrons with the electron spin fluctuations antisymmet- 
ric with respect to the layers. Such fluctuations correspond 
to "acoustical" gapless spin-wave branch in the antiferromag- 
netic insulating phase. This is natural, since other fluctuations 
immediately decay into quasiparticles. On the other hand, 
such fluctuations are rather long ranged, and therefore they 
define, actually, the true "coherent" gap. It is no surprise, 
therefore, that with underdoping the position of the neutron 
maximum varies proportional to Tc and not to the gap in the 
quasiparticle spectrum measured in tunneling experiments. 



The chemical potential /ii is calculated with respect to the 
bottom of the band. The term with ar is negative below Tc 
and it keeps '^a.b attached to a certain vicinity, whereas the 
term with arj corresponds to hopping between them. With in- 
creasing \t\ the motion becomes more one-dimensional, and 
eventually the orbital motion is unable to destroy supercon- 
ductivity, i.e., Hc2- For sufficiently small values of h this can 
happen even close to Tc (see Fig. 11). In the immediate vicin- 
ity of Tc 




For sufficiently high fields another mechanism becomes im- 
portant, namely the paramagnetic destruction of Cooper pairs. 
Unfortunately, the account of this effect leads to the restora- 
tion of a negative curvature at lower temperatures and larger 
fields. This does not fit the experimental data (see [^^). 

There was an interesting suggestion [|6|| that in the cases 
of reduced dimensionality (films, wires in longitudinal field), 
where the orbital destruction becomes impossible, the param- 
agnetic effect leads to a sequence of phase transitions between 
different inhomogeneous LOFF (Larkin-Ovchinnikov-Fulde- 
Ferrell) phases with increasing critical fields. This is a possi- 
ble origin of the positive curvature at very low temperatures. 



VIII. UPPER CRITICAL FIELD 
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The last interesting consequence of the model, which I will 
describe, is the positive curvature of the upper critical mag- 
netic field, as function of temperature [ p3[ | which was ob- 
served in many experiments, (e.g. [p4|], see Fig. 10). The usual 
calculation gives a negative curvature. The reason for such 
a change is contained in the form of the spectrum. As it 
was mentioned before, the most important are the vicinities of 
the singular regions, where the motion of electrons becomes 
one-dimensional. There is, however, a hopping between the 
"a" and "b" vicinities, and hence a two-dimensional motion. 
Therefore, close to Tc the field Hc2 is defined by the usual 
orbital motion. With departing from Tc the situation changes. 
This can be seen from the Ginsburg-Landau type equations 
derived in the vicinity of Tc- They correspond to the mini- 
mum of the free energy 



dx 



. d 2e 

dx c 
d 2e 



where 



2(7rT,)^ 
7C(3)mi' 



T-Tc 



Tc 



-,v = 



IA2I 



(Ai + |A2|)= 



(21) 
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FIG. 4. ARPES result for the dependence of quasiparticle energy 
on components of quasimomentum in the vicinity of the extended 
saddle point singularity. 




FIG. 1. Experimental results for the isotope shift (Oie — > Ois) 
in partially substituted YBCO {Y — > Pr, and Ba -» La). The 
exponent a in Tc oc M°' is presented as function of Tc together with 
the linear regression. 



ens: * — ;: \<i-^ 
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P ^5 36 ^5 M 

FIG. 2. Compilation of tunneling a^F^uo) for NCCO (solid line) 
and the phonon density of states F{lo) generated from the disper- 
sion curves obtained from the inelastic neutron scattering from sin- 
gle-crystal NCCO 




FIG. 5. Fermi surface with singular regions having a high density 
of states. A short range interaction connects any point with the sin- 
gular region, whereas a long range interaction connects only points 
within this region. 




FIG. 6. Plot of A/Ai according to eq. |. 






/ 
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FIG. 3. The same, as in Fig. 2, for BSCCO. The lower curve is 
F(lo) from inelastic neutron scattering 



FIG. 7. Schematic presentation of the superconducting bands in 
the vicinity of singular points. The abscissa is p — a:: for the a- vicinity, 
and py for the b- vicinity. The notations 1,1', 2,2' correspond to ex- 
tremal points with the maximal density of states. 
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FIG. 8. Plot of the imaginary part of the electronic spin suscepti- 
bility, which is the energy dependent factor in the neutron inelastic 
scattering cross section Xm x (m-effective mass, d-period along), 
as function oi 8 = (SI/ A — 2)^'^^ where Q, is the energy loss, in the 
case 7 = /xi/A = 0.1 (/ii is the Fermi energy with respect to the 
bottom of the ID band). Here the point f2 = 2A is the threshold, 
and below this poind x" = 



Y3a CuO 

I I HI 



FIG. 9. Experimental data from for the same quantity in ar- 
bitrary units as function of energy loss for optimally doped YBCO. 
The initial rise corresponds to the smeared out threshold. 




FIG. 10. Experimental results 
dence of the upper critical field 



for the temperature depen- 
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